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Abstract 


Recent  progress  is  reported  in  development  of  ab  initio  computational  methods  for  the 
electronic  structures  of  molecules  employing  atomic  spectral-product  representations.  In 
this  approach,  the  physically  significant  many-electron  eigenstates  of  constituent  atoms 
provide  a  computational  basis  for  descriptions  of  the  electronic  structure  of  matter  as  an 
alternative  to  more  commonly  employed  atomic-  or  molecular-orbital-based  representa¬ 
tions.  The  Hamiltonian  matrix  is  seen  to  be  a  sum  over  atomic  energies  and  a  pairwise 
sum  over  Coulombic  interaction  terms  which  depend  only  on  the  separations  of  the  in¬ 
dividual  atomic  pairs.  Overall  electron  antisymmetry  can  be  enforced  in  the  approach 
by  unitary  transformation  when  appropriate,  rather  than  as  a  possibly  encumbering  or 
unnecessary  global  constraint.  The  matrix  representative  of  the  antisymmetrizer  in  the 
spectral-product  basis,  which  is  equivalent  to  the  metric  matrix  of  the  corresponding  ex¬ 
plicitly  antisymmetric  basis,  provides  the  required  transformation  to  antisymmetric  or 
linearly-independent  states  after  Hamiltonian  evaluation.  Particular  attention  is  focused 
in  the  present  report  on  properties  of  the  metric  matrix  and  on  the  spectral  composi¬ 
tions  of  molecular  eigenstates  described  in  spectral-product  representations.  Illustrative 
calculations  are  reported  for  simple  but  protot-ypically  important  diatomic  (H2,  CH)  and 
triatomic  (H3,  C 1 1 2 )  molecules  employing  recently  devised  algorithms  and  computer  codes. 
This  particular  implementation  of  the  approach  combines  valence-bond  constructions  of 
standard  tableau  functions  and  matrices  with  transformations  to  atomic  eigenstate  product 
representations,  employing  Slater-orbital-based  one-  and  two-electron  integral  evaluations. 
The  calculated  metric  matrices  and  corresponding  potential  energy  surfaces  obtained  in 
this  way  elucidate  a  number  of  aspects  of  the  spectral-product  development,  including  the 
nature  of  closure  in  the  spectral-product  representation,  the  general  redundancy  of  its  ex¬ 
plicitly  antisymmetrized  form,  the  convergence  of  the  two  apparently  disparate  descriptions 
to  a  common  invariant  subspace,  and  the  manner  in  which  a  chemical  bonding  descriptor 
is  provided  by  the  atomic  spectral  compositions  of  molecular  eigenstates  obtained  in  the 
representation.  Concluding  remarks  describe  briefly  additional  studies  in  progress  and  the 
prognosis  for  performing  spectral-product  calculations  more  generally  and  efficiently. 
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I.  Introduction 


The  electronic  eigenstates  of  atoms  provide  a  well-known  formal  basis  for  simple  product 
representations  of  molecular  electronic  wave  functions  and  interaction  energies  in  long- 
range  separation  limits.1  More  generally,  ab  initio  investigations  of  the  electronic  structures 
and  associated  potential  energy  surfaces  of  molecules  tend  to  employ  explicitly  antisym¬ 
metric  representations  based  on  early  molecular-orbital  or  valence-bond  descriptions  of  the 
relevant  electronic  degrees  of  freedom.2-5  In  these  approaches,  the  presence  of  the  compo¬ 
nent  atomic  constituents  is  not  made  explicit,  but  is  evident  largely  in  the  choices  of  atomic 
orbitals  employed  in  the  molecular  representational  basis  sets.  It  has  apparently  proved 
difficult  to  base  a  general  ab  initio  computational  approach  to  the  electronic  structures 
of  matter  directly  on  an  atomic  representation,6  in  spite  of  the  fact  that  atoms  clearly 
comprise  molecules  and  other  forms  of  matter. 

A  recent  series  of  reports  emphasizes  both  the  difficulties  and  potential  advantages  of 
adopting  an  atomic  eigenstate-product  representation  in  ab-initio  calculations  of  molecular 
electronic  structures.7-11  This  program  of  study,  the  motivations  and  purposes  of  which 
were  described  earlier,12  is  predicated  largely  on  the  demonstrable  closure  of  the  outer 
spectral  product  of  complete  constituent  atomic  eigenstates  for  representations  of  molecular 
electronic  states  in  the  absence  of  explicit  (term-by-term)  enforcement  of  aggregate  electron 
antisymmetry.  Although  both  conceptual  and  computational  barriers  must  be  overcome  in 
such  an  approach,  certain  potential  advantages  follow  from  its  adoption.  These  include  the 
possibility  of  performing  accurate  atomic  and  atomic-interaction  calculations  once  and  for 
all  and  retaining  such  information  for  repeated  molecular  applications,  employing  entirely 
analytical  methods  in  calculating  the  angular  dependences  of  atomic  interactions,  and  the 
opportunity  to  enforce  molecular  antisymmetry  when  required,  rather  than  as  a  possibly 
encumbering  overall  global  constraint,  to  mention  some  important  examples.'-12 

In  the  present  contribution,  additional  theoretical  aspects  of  the  atomic  spectral-product 
approach  to  molecular  electronic  structure  are  reported,  the  essential  features  of  a  computer 
code  suite  for  performing  such  calculations  are  described,  and  selected  computational  appli¬ 
cations  to  prototypically  important  diatomic  (H2,  CH)  and  triatomic  (H3,  CH2)  molecules 
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are  presented.  Calculations  of  metric  matrices  in  these  cases,  and  of  corresponding  po¬ 
tential  energy  curves  and  surfaces,  in  spectral-product  representations  illustrate  the  capa¬ 
bilities  of  the  valence-bond-  and  Slater-orbital-based  algorithms  and  codes  devised.  The 
metric  matrices,  which  provide  a  connecting  bridge  between  spectral-product  representa¬ 
tions  in  the  absence  or  presence  of  term-by-term  antisymmetry,  clarify  a  number  of  aspects 
of  the  development,  including  the  nature  of  closure  in  the  spectral-product  representation 
and  the  general  redundancy  of  its  explicitly  antisymmetrized  form.  The  calculations  re¬ 
ported  for  H2  and  H3  characterize  the  spectrum  of  the  antisymmetrizer  as  that  of  a  compact 
operator  in  the  spectral-product  representation,  and  demonstrate  specific  redundancies  be¬ 
tween  commonly  employed  charge-transfer  and  one-electron  covalent  excitations,  whereas 
the  potential  energy  curves  for  ground  and  excited  states  demonstrate  the  advantages 
of  the  Slater-orbital-based  integral  evaluations  employed.  The  spectral  compositions  of 
the  atomic-product  representations  of  molecular  eigenfunctions  are  seen  to  provide  highly 
satisfying  descriptors  of  chemical  bonding  in  CH  and  CH2  molecules,  which  have  played 
important  historical  roles  in  the  development  of  accurate  ab  initio  quantum  methods  for 
electronic  structure  and  spectra.13,14 

The  theory  is  described  in  Section  II,  where  the  atomic  spectral  representations  of  aggre¬ 
gate  electronic  degrees  of  freedom  are  defined,  the  forms  of  metric  and  Hamiltonian  ma¬ 
trices  with  and  without  explicit  antisymmetrization  are  reported,  and  the  unitary  trans¬ 
formation  formalism  for  isolation  of  the  physical  and  linearly-independent  Hamiltonian 
subspaces  are  given  and  their  equivalence  established.  Calculations  employing  the  recently 
devised  algorithms  and  computer  codes,  reported  in  Section  III,  include  metric  matrices 
and  energy  eigenvalues  for  the  aforementioned  H2,  CH,  H3,  and  CH2  molecules.  Discussion 
and  concluding  remarks  presented  in  Section  IV  describe  additional  studies  in  progress  and 
strategies  employed  in  performing  atomic  spectral-product  calculations  more  generally  and 
efficiently. 
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II.  Theoretical  Background  and  Developments 

Selected  background  information  is  provided  in  this  Section  in  summary  of  previously 
devised  aspects  of  the  spectral-product  approach  to  molecular  electronic  structure,7-12 
and  new  theoretical  developments  and  clarifications  are  reported.  The  atomic  spectral- 
product  representations  of  molecular  electronic  states  of  interest  here  are  described  in 
Section  A,  construction  of  Hamiltonian  matrices  in  these  representations  are  described  in 
Section  B,  and  the  role  of  the  aggregate  metric  matrix  in  isolation  of  physically  significant 
Schrodinger  eigenstates  is  described  in  Section  C. 


A.  Atomic  Spectral-Product  Representations. 

The  adiabatic  (Born-Oppenheimer)  electronic  states  obtained  from  solution  of  the  noil- 
relativistic  Schrodinger  equation 

H(r  :  R)^f(r  :  R)  —  \I>(r  :  R)  ■  E(J2)  (1) 


provide  useful  first  approximations  to  the  chemical  and  physical  attributes  of  molecules  and 
other  atomic  aggregates.15  Here,  the  row  vector  vE'(r  :  R)  contains  the  desired  eigenstates, 
the  diagonal  matrix  E (R)  contains  the  corresponding  energy  eigenvalues,  r  refers  to  the 
space  and  spin  coordinates  of  the  rp  electrons  of  the  system,  and  R  designates  the  fixed 
positions  (Ri,  R2, . . .  Rn)  of  the  N  atomic  nuclei. 


The  non-relativistic  Hamiltonian  operator  in  Eq.  (1), 

-2  N  nt  ry  _  2  nt  —  lnt 

X 

2m 


nt 


6(r:K>  =  'Z-£;*l-t£¥-+'E  £ 


Tkc 


Tkk’ 


N- 1  N 

+  £  £ 


ZaZpe2 


R, 


a/3 


(2) 


fc= 1  a=l  k=l  n',x  fc= 1  k'=k-\-l  a=l  /3=a+l 

contains  the  familiar  kinetic,  electron-nuclear  attraction,  and  electron  and  nuclear  repulsion 
operators  and  is  totally  symmetric  in  all  electron  coordinates. 16  Consequently,  the  solutions 
of  Eq.  (1)  can  be  classified  according  to  the  irreducible  representations  of  the  symmetric 
group  Snt,  with  the  totally  antisymmetric  solutions  providing  the  physically  significant 
Schrodinger  eigenstates.1' 

form7-10 


The  atomic  spectral-product  representation  is  written  in  the  outer-product 
$(r  :  R)  -  j$(1)(i)  ®  $(2)(2)  ®  •  •  •  $(A°(n)  , 


(3) 
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where  the  row  vectors  3>(a)(z)  contain  the  orthonormal  antisymmetric  eigenstates  of  the 
constituent  atoms  (a  =  1  to  N )  specihed  by  the  quantum  numbers  (E.  L,  Ml ,  S.  Ms ,  P)a,18 
with  i  referring  to  the  coordinates  of  the  na  electrons  arbitrarily  assigned  to  the  atom  a 
measured  relative  to  the  atomic  centers  Ra.  The  subscript  “O”  in  Eq.  (3)  indicates  the 
adoption  of  a  particular  ordering  convention  for  the  vector  sequence  of  spectral-product 
functions  employed.  More  detailed  descriptions  of  these  and  of  the  other  notational  con¬ 
ventions  of  Eqs.  (1)  to  (3)  are  reported  elsewhere.'-12 


Although  the  representation  of  Eq.  (3)  is  not  term-by-term  antisymmetric  in  all  electron 
coordinates,  it  is,  in  fact,  complete  in  the  limit  of  spectral  closure  for  representations  of 
totally  antisymmetric  aggregate  eigenstates.1,19  Accordingly,  Eq.  (3)  is  in  this  connection 
formally  equivalent  to  the  more  familiar  explicitly  (term-by-term)  antisymmetric  form6 


:  R)  =  pA&ir  ■  R ) 


=  PA  j#(1)(i)  <g>  $(2)(2)  ®  •  •  •  $(Ar)(n)|o  , 


(4) 


where 


nt\ 

Pa  =  {nt\)~1/2{n i!n2!  ■  •  -nNl)~1/2  ^(-l)5pPp  (5) 

p=i 

is  the  aggregate  antisymmetrizer,17  normalized  to  insure  that  (PA  3?(r*  :  R)  |  P_,i  ^(r  :  R  j)  — > 
I  in  the  limit  (R  — >  oo),  in  accordance  with  the  orthonormality  of  the  square-integrable 
antisymmetric  atomic  eigenstates  in  the  row  vectors  Alternatively,  it  is  possible 

to  enforce  aggregate  electron  antisymmetry  employing  a  coset  decomposition  of  the  anti¬ 
symmetrizer  to  isolate  those  electron  permutations  not  completely  internal  to  the  already 
antisymmetric  atomic  eigenstates.6,20 


An  important  difference  between  the  foregoing  two  representations  is  found  in  their  re¬ 
spective  closure  relations.  In  the  case  of  the  orthonormal  representation  of  Eq.  (3),  the 
closure  expression  is  <I?(r  :  R)  ■  $(r'  :  Ry  =  S(r  —  r'  :  R),  where  S(r  —  r'  :  R)  applies  to 
all  irreducible  representations  of  Snt  spanned  by  the  basis,  whereas  the  closure  expression 
is  :  R )  ’  *&a{v  :  Ry  —  Q^A(r  —  r'  :  R)  for  the  non-orthogonal  explicitly  antisym¬ 

metric  representation  of  Eq.  (4),  where  Q  —  nt\l(n\\ny.  ■  ■  -njv!)  is  the  redundancy  of  the 
explicitly  antisymmetric  basis  of  Eq.  (4)  in  the  closure  limit.10 
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B.  Hamiltonian  Matrices  in  the  Spectral-Product  Representations. 


It  is  convenient  to  rewrite  the  Hamiltonian  operator  of  Eq.  (2)  in  a  form  suggested  by 
the  arbitrary  assignments  of  electrons  to  atoms  made  in  Eq.  (3)  in  constructing  spectral- 
product  matrix  representatives.  Specifically,  the  Hamiltonian  operator  of  Eq.  (2)  is  written 
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includes  the  electronic  terms  arbitrarily  associated  with  the  atom  a,  and 
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13 


(6) 

(7) 

(8) 


includes  the  Coulombic  interaction  terms  associated  with  the  pair  of  atoms  (a,/3). 


Employing  the  basis  sets  of  Eqs.  (3)  and  (4)  in  variational  solutions  of  Eq.  (1)  gives  the 
matrix  Schrodinger  equations 


H (R)  •  Uh(-R)  -  Uh(-R)  ■  E (R),  (9) 

and 

H a(R)  ■  VH (R)  =  S (R)  •  Vh(-R)  •  E(.R),  (10) 

where 

N  N-l  N 

H(R)  =  {*(r:R)\H(r:R)\*(r:R))  =  Y,  H(a)  +  ^  £  V^(Ra/3),  (11) 

a=l  a=l  /3=a+ 1 

and 

N  N-l  N 

HA(R)  =  (^A(r:R)\H(r:R)\^A(r:R))  =  J2  H<“>(fl)+^  £  V^(R), 

a= 1  a=l  /3=a+l 

(12) 

are  the  respective  Hamiltonian  matrices. 
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In  Eq.  (9),  the  unitary  solution  matrix  Uh(-R)  diagonalizes  the  Hamiltonian  matrix  of 
Eq.  (11)  constructed  in  the  orthonormal  basis  of  Eq.  (3).  By  contrast,  the  matrix  Vh(-R) 
in  Eq.  (10)  is  generally  not  unitary  since  the  metrically-defined  Hamiltonian  matrix  of  Eq. 
(12)  is  constructed  in  the  lion-orthogonal  basis  of  Eq.  (4).  Rather,  Vh(H)  satisfies  the 
condition  Vh(-R)'  •  S (R)  ■  'Vn(R)  —  I,15  with  the  metric  matrix  S (R)  here  and  in  Eq. 
(10)  given  by  the  expression 

S (R)  =  (PA®(r  :  R)\PA<!>(r  :  R))  =  Q1/2^  :  R)\PA\<f>(r  :  R)),  (13) 

where  the  factor  Q  =  rp!/(ni!n2!  ■  ■  •  njv!)  is  the  aforementioned  redundancy  of  the  explicitly 
antisymmetric  basis  of  Eq.  (4)  in  the  closure  limit.10  Evidently,  S(R)  has  an  interpretation 
as  the  matrix  representative  of  the  antisymmetrizer  of  Eq.  (5)  in  the  basis  of  Eq.  (3), 
separate  from  its  role  as  the  metric  matrix  of  the  noil-orthogonal  basis  of  Eq.  (4). 

Although  the  Hamiltonian  matrices  of  Eqs.  (11)  and  (12)  are  similar  in  form,  the  individual 
atomic  and  atomic-pairwise  interaction  terms  shown  there  differ  significantly  in  the  two 
representations.  Specifically,  the  atomic  terms  in  Eq.  (11)  obtained  from  Eq.  (7)  are 
independent  of  atomic  position,  and  the  Coulombic  interaction  terms  obtained  from  Eq. 
(8)  depend  only  on  the  vector  separation  R^p  of  an  individual  atomic  pair  (a,/3).  By 
contrast,  the  terms  in  Eq.  (12),  which  must  be  evaluated  using  the  Hamiltonian  of  Eq. 
(2),  formally  depend  on  the  positions  R  of  all  the  atoms  in  the  aggregate  due  to  the  overall 
electron  antisymmetry  of  the  basis  of  Eq.  (4).  That  is,  because  the  electronic  coordinates 
in  the  representation  of  Eq.  (4)  appear  at  all  the  atomic  centers,  the  individual  Coulombic 
interactions  in  the  Hamiltonian  of  Eq.  (2)  can  be  alternatively  inter-  or  intra-atomic  terms, 
and  the  partitioning  of  Eqs.  (6)  to  (8),  which  constitutes  an  assignment  of  particular 
electrons  to  specific  atomic  centers,  can  not  be  employed  in  this  case. 

The  simpler  forms  of  the  terms  appearing  in  Eq.  (11)  relative  to  those  in  Eq.  (12)  are 
largely  consequences  of  the  orthonormality  of  the  spectral-product  basis  employed  in  the 
absence  of  prior  enforcement  of  overall  aggregate  electron  antisymmetry,  the  use  of  atomic 
eigenstates  in  the  representation,  and  the  atomic  pairwise-additive  nature  of  the  interaction 
terms  in  the  Hamiltonian  operator  of  Eq.  (6).  By  contrast,  the  atomic-pairwise  nature 
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of  the  interaction  terms  in  the  Hamiltonian  operator  can  not  be  employed  to  advantage 
in  conjunction  with  use  of  the  explicitly  antisymmetric  representation  of  Eq.  (4).  It  spite 
of  the  apparent  differences  in  the  two  representations,  both  can  give  the  same  physical 
Schrodinger  eigenstates  \J/p(r  :  R)  =  {$(r  :  R)  ■  Uu{R)}P  —  {3?a(»"  :  R)  •  Vh(-R)}p  and 
energies  {E(i?)}pp  under  appropriate  conditions. 

C.  Isolation  of  Physical  Eigenstates  in  the  Spectral-Product  Representations. 

It  has  been  shown  that  the  basis  of  Eq.  (3)  spans  the  totally  antisymmetric  representation 
of  the  aggregate  rq-electron  symmetric  group  Sri  f  once  and  only  once  in  the  limit  of  closure, 
although  other  unphysical  non-totally-antisymmetric  (non-Pauli)  representations  of  Snt 
are  also  spanned  by  the  basis.10  The  non-Pauli  solutions  are  accordingly  included  in  the 
columns  of  the  unitary  transformation  matrix  Uh(-R)  of  Eq.  (9).  By  contrast,  although 
the  solutions  of  Eq.  (10)  are  totally  antisymmetric  by  construction,  they  generally  become 
unstable  in  large  representations,  consequent  of  the  linear  dependence  arising  from  the 
Q-fold  redundancy  of  the  basis  of  Eq.  (4)  in  the  closure  limit.  These  two  assertions  are 
closely  related,  and  suggest  the  possible  equivalence  of  results  obtained  employing  either 
“prior”  or  “post”  antisymmetrization  in  extracting  physically  acceptable  solutions  from 
the  molecular  Schrodinger  equation  in  the  two  representations  of  Eqs.  (3)  and  (4). 

The  desired  Pauli  and  linearly-independent  solutions  can  be  obtained  by  transforming 
Eqs.  (9)  and  (10)  to  representations  which  isolate  totally  antisymmetric  and  linearly- 
independent  subspaces,  respectively,  and  can  be  shown  to  be  equivalent  in  the  limit  of 
spectral  closure.  The  unitary  transformation  matrix  Us(-R)  required  to  accomplish  this 
partitioning  of  Eqs.  (9)  and  (10)  is  obtained  from  diagonalization 

Us(fi)*.S(«).UsW  =  Sd(fl)=(W«)>-  o)  <14> 

of  the  Hermitian  metric  matrix  of  Eq.  (13).  Here,  the  non- negative  eigenvalue  matrix 
s d(R)  is  partitioned  into  an  upper  diagonal  {sf/(R)}pp  containing  the  largest  eigenvalues 
ordered  in  decreasing  value,  which  tend  to  the  indicated  common  redundancy  factor  Q  in 
the  closure  limit,  and  a  lower  diagonal  {sd(-R)}WM  of  eigenvalues  which  tend  to  zero  in  this 
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limit,  so  ordered  by  appropriate  arrangement  of  the  columns  of  the  transformation  matrix 
Us(-R)  into  physical  ( p )  and  unphysical  ( u )  blocks. 

The  diagonalization  of  Eq.  (14)  is  seen  to  correspond  to  construction  of  eigenstates  of  the 
Hermitian  antisymmetrizer  in  the  spectral-product  basis  of  Eq.  (3)  in  the  form 

#s(r-  :  R)  =  $(r  :  R )  •  US(R)  -  {{$s(r  :  R)}p,  :  i2)}w},  (15) 

where  { (n  :  R)}p  contains  totally  antisymmetric  states  corresponding  to  the  non-zero 
eigenvalues  of  the  antisymmetrizer,  and  { <E>  s  ( r  ■  R)}u  contains  non-Pauli  states  corre¬ 
sponding  to  the  zero  eigenvalues  of  the  antisymmetrizer  in  the  closure  limit.  By  contrast, 
in  the  explicitly,  or  term-by-term,  antisymmetric  basis  of  Eq.  (4),  the  states  generated  by 
the  transformation  of  Eq.  (14)  are  written  in  the  canonically  orthogonalized  form 

$s (r  :  R)  =  #a(t*  :  R)  •  U S(R)  ■  s^R)-1/2  -  {{$s(r-  :  R)}P,  {*s (r  :  R)}u  -  0},  (16) 

which  are  correctly  normalized  linearly-independent  combinations  of  the  original  basis  and 
un-normalizable  null  states  associated  with  linearly-dependent  combinations  of  the  basis 
of  Eq.  (4). 21  In  the  case  of  degenerate  eigenvalues,  the  eigenstates  of  Eqs.  (15)  and  (16) 
are  generally  arbitrary  but  are  nevertheless  separated  into  invariant  subspaces.  Of  course, 
in  the  large-separation  limit  S(R  — >  oc)  — »  I,  in  which  case  the  eigenvalues  of  the  metric 
matrix  are  all  unity,  both  representations  of  Eqs.  (3)  and  (4)  are  orthonormal  and  the 
states  of  Eqs.  (15)  and  (16)  are  not  required. 

Equations  (14)  and  (15)  indicate  that  the  spectrum  of  PA  acting  in  the  domain  of  the 
spectral-product  basis  of  Eq.  (3)  is  that  of  a  compact  Hermitian  operator,22  with  zero 
providing  a  lower  limiting  point  of  accumulation  of  the  spectrum  and  the  associated  non- 
Pauli  states  {<f»s(^  :  R)}u  of  Eq.  (15)  corresponding  to  its  null  space  Pa{^s(^  :  R)}u  — > 
0.  The  states  {3?s(r  :  R)}p  of  Eq.  (15)  are  associated  with  the  upper  limiting  point, 
PA{®s (r  :  R)}P  -  Q1/2{<Mr-  :  R)}p,  and  correspond  to  the  physically  significant  totally 
antisymmetric  spectrum  of  the  antisymmetrizer.  The  linearly-independent  combinations 
of  the  term-by-term  antisymmetric  basis  functions  of  Eq.  (4)  given  by  Eq.  (16)  are 
identical  to  the  totally  antisymmetric  states  of  Eq.  (15)  in  the  closure  limit,10  whereas  the 
linearly-dependent  states  of  Eq.  (16)  are  seen  to  be  eliminated  by  the  development.21 
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The  physically  significant  blocks  of  the  Hamiltonian  matrices  of  Eqs.  (11)  and  (12)  are 
obtained  from  the  matrix  Us(-R)  in  the  upper  left-hand  forms 

{H(i?)}pp  =  {Us(JR)t  •  H(R)  ■  Us (R)}pp  (17) 

{Ua(R)}pp  =  {s d(R)-1/2}pp  •  {Us(i?)f  •  H a(R)  •  Us (R)}PP  •  {s d(R)-1/2}pp,  (18) 

where  now  all  terms  in  both  Hamiltonian  matrices  depend  formally  on  the  positions  R  of 
all  the  atoms  in  the  aggregate.  The  unitary  transformation  matrix  Us(-R)  of  Eq.  (14) 
obtained  from  the  metric  matrix  of  Eq.  (13)  is  seen  to  be  responsible  for  incorporating  the 
non-local  effects  of  inter-atomic  aggregate  electron  permutation  symmetry  in  the  Hamil¬ 
tonian  matrix  of  Eq.  (11)  by  virtue  of  Eq.  (17).  Overall  electron  antisymmetry  can  be 
enforced  as  required  in  this  approach  in  accordance  with  the  spatial  separations  of  the  in¬ 
dividual  atoms  and  the  perceived  strengths  of  their  interactions  in  the  Hamiltonian  of  Eq. 
(17),  rather  than  as  a  possibly  encumbering  overall  global  constraint,  as  employed  in  the 
Hamiltonian  matrix  in  Eq.  (18).  Although  the  role  of  the  transformation  matrix  Us(-R)  in 
the  explicitly  antisymmetric  representation  of  Eq.  (4)  is  to  isolate  the  linearly-independent 
subspace  present  therein,  it  would  seem  useful  to  also  employ  approximation  methods  in 
its  evaluation  in  certain  cases,  although  the  more  significant  advantage  of  this  approach  is 
expected  to  be  found  in  connection  with  Eq.  (17). 

Finally,  identical  physically  significant  Schrodinger  eigenstates  are  obtained  from  the  imme¬ 
diately  foregoing  Hamiltonian  matrices  in  the  form  T' p(r  :  R)  =  {$(r  :  R)}{>  •  {Uh(-R)}pp, 
where  {$(r  :  R)}p  is  obtained  from  either  Eq.  (15)  or  (16)  and  {Uh(-R)}pp  is  obtained 
from  either  of  the  Hamiltonian  matrices  of  Eqs.  (17)  or  (18).  Although  the  requirements 
of  computational  implementations  of  the  two  approaches  differ  significantly,  in  view  of  the 
“post”  and  “prior”  incorporation  of  aggregate  electron  antisymmetry  in  the  two  Hamilto¬ 
nian  matrices,  the  developments  are  clearly  united  through  use  of  the  metric  matrix  of  Eq. 
(13)  in  isolating  a  common  invariant  physical  subspace  in  which  to  construct  Schrodinger 
eigensolutions. 
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III.  Computational  Implementations  and  Applications 


Algorithms  and  computer  codes  have  been  devised  to  evaluate  the  various  expressions 
given  in  the  preceding  Sections.  In  brief,  a  valence-bond  code  suite  (CRUNCH23)  has  been 
modified  to  accommodate  molecular  integrals  evaluated  in  real  Slater  orbitals  (SMILES24), 
and  transformations  devised  from  the  standard  tableau  descriptions  of  configurational  state 
functions  and  Hamiltonian  and  metric  matrix  elements  of  the  valence-bond  formalism  to 
the  atomic  spectral-product  expressions.  Computational  results  obtained  in  this  way  are 
reported  in  Section  A  for  H2  and  H3  molecules,  in  Section  B  for  diatomic  CH,  and  in 
Section  C  for  triatomic  CH2.  These  calculations  are  illustrative  of  selected  aspects  of  the 
approach,  and  of  the  capabilities  of  the  computational  algorithms  and  codes  devised,  and 
are  not  meant  to  provide  exhaustive  computational  descriptions  of  the  quantities  reported. 

A.  Diatomic  and  Triatomic  Hydrogen  Molecules. 

Although  the  metric  matrix  is  familiar  in  connection  with  canonical  transformations  of 
one-electron  orbital  basis  sets,21  its  attributes  and  importance  in  connection  with  many- 
electron  basis  states  are  perhaps  less  familiar,  particularly  in  the  context  of  the  spectral- 
product  representations  described  here.  Of  particular  interest  is  the  function  of  the  metric 
matrix  in  providing  a  measure  of  the  spectral  closure  of  the  representation  of  Eq.  (3)  in 
the  absence  of  explicit  antisymmetry  and  of  the  redundancy  of  the  representation  of  Eq. 

(4) ,  particularly  in  connection  with  commonly  employed  ionic  terms  involving  the  transfer 
of  one  or  more  electrons  from  one  atom  to  others.  Illustrative  molecular  calculations 
involving  hydrogen  atoms  help  to  clarify  these  issues,  and  demonstrate  the  capabilities  of 
the  algorithms  and  codes  devised. 

In  Figure  1  are  shown  the  eigenvalues  of  the  metric  matrix  of  Eq.  (13)  for  symmetric 
collinear  H3  as  a  function  of  atomic  separation  (Rab  =  Rbc),  constructed  in  (a)  2slp 
and  (b)  3s2p  atomic  Slater  basis  sets  employing  hydrogenic  exponents.  These  atomic 
orbitals  give  representations  <jT-H3)(r  :  R )  =  |$(H“)(i)  ®  #(^>(2)  ®  3?(Hc)(3)}0  which 
include  375  and  2187  terms,  respectively.  The  eigenvalues  of  the  antisymmetrizer  of  Eq. 

(5)  in  these  representations  are  seen  to  appear  only  in  the  allowable  interval  (0,(5)  where 
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Q  =  3!/(  1!  1!  l!)=6,  and  to  fill  the  interval  for  smaller  values  of  atomic  separation.  The  great 
many  values  near  0  in  Panel  (b)  are  in  accordance  with  the  lower  end  of  the  spectrum  of  Pa 
providing  a  point  of  accumulation.  The  eigenstates  corresponding  to  these  approximately 
zero  eigenvalues  refer  to  non-totally  antisymmetric  states  in  the  product  basis  [Eqs.  (3)  and 
(15)],  whereas  in  the  prior  antisymmetrized  basis  they  refer  to  linearly-dependent  states 
[Eqs.  (4)  and  (16)].  Conversely,  the  states  of  Eqs.  (15)  and  (16)  corresponding  to  the 
largest  eigenvalues  of  Figure  1  refer  to  approximately  totally  antisymmetric  and  linearly- 
independent  states,  respectively.  Note  that  the  value  Q= 6  implies  the  antisymmetric  basis 
of  Eq.  (4)  is  six-fold  redundant  in  the  limit  of  closure  in  this  case,  which  redundancy  is 
removed  by  the  canonical  orthogonalization  procedure  of  Eq.  (16).  As  the  number  of  basis 
functions  and  their  spatial  extent  are  increased,  the  eigenvalues  of  Figure  1  converge  to 
the  upper  and  lower  limiting  points  of  the  spectrum  for  finite  values  of  atomic  separation, 
whereas  they  approach  unity  in  the  limit  of  large  separation  in  any  square-integrable 
representation. 

Although  large  orbital  basis  sets  are  generally  required  to  achieve  convergence  in  spectral- 
product  representation,  metric  matrices  constructed  even  in  small  basis  sets  provide  useful 
information.  In  Figure  2  are  shown  constant- value  contours  for  the  two  doublet  eigenvalues 
of  the  metric  matrix  for  T-shaped  H3  as  functions  of  width  (w)  and  height  ( h )  obtained  in  a 
minimal  Is3  spectral- product  representation.  In  panel  (a),  the  contours  are  monotonically 
increasing  from  that  labeled  1.0,  whereas  in  panel  (b)  they  are  monotonically  decreasing 
from  that  labeled  0.9.  The  cusps  in  the  contours  of  both  panels  indicate  the  presence 
of  a  seam  of  intersection  in  the  two  surfaces  along  the  w  —  1.155  h  line  corresponding 
to  1)?,},  symmetry.  This  high-symmetry  seam  corresponds  to  the  better-known  seam  of 
intersection  of  the  two  lowest-lying  doublet  energy  surfaces  in  H3,25  identified  here  entirely 
on  basis  of  the  eigenvalues  of  the  metric  matrix  in  the  absence  of  energy  calculations.  Such 
calculations  of  metric  matrices,  which  entail  evaluations  only  of  overlap  matrix  elements, 
can  provide  useful  information  in  other  molecular  cases  more  generally. 

The  spectral  product  representation  of  Eq.  (3)  is  complete  employing  only  the  indicated 
neutral  atomic  states,  and  the  associated  antisymmetrized  basis  of  Eq.  (4)  is  correspond- 
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ingiy  Q-fold  redundant  in  the  closure  limit.  It  is  instructive,  however,  to  examine  the 
effects  of  adding  commonly  employed  ionic  or  charge-transfer  configurations  to  the  ba¬ 
sis.  In  Figure  3  are  shown  the  singlet  and  triplet  eigenvalues  of  the  metric  matrix  of  Eq. 
(13)  for  diatomic  hydrogen  constructed  in  a  7s5p3d2flg  Slater  basis  set  employing  fixed 
hydrogen  Is,  2 p,  3d ,  4/,  5g  exponents  in  Sturmian  sequences.26  The  eigenvalues  shown  cor¬ 
respond  to  single-excitation  configuration-interaction  calculations,  including  both  covalent 
(H-H  -  119  terms)  and  ionic  (H+H_  and  H+H_  -  154  terms  )  structures.  Although  the 
redundancy  factor  is  Q  —  2!/ (1!1!)=2  for  the  spectral-product  or  covalent  representations 
of  Eqs.  (3)  and  (4)  in  this  case,  the  presence  of  the  additional  charge-transfer  terms  gives 
rise  to  eigenvalues  which  appear  in  the  larger  interval  (0,4).  This  is  a  consequence  of  the 
four-fold  redundancy  of  the  combined  covalent  and  ionic  basis  employed  in  the  calculations 
of  Figure  3  in  the  closure  limit.  Specifically,  the  band  of  near  zero  eigenvalues  in  Figure  3 
arises  from  the  near  linear-dependence  of  charge-transfer  terms  with  the  diffuse  covalent 
excitations  -  this  band  is  absent  in  the  eigenvalue  spectrum  (0,2)  obtained  with  covalent 
terms  alone.9,10  Moreover,  the  bands  of  eigenvalues  approaching  2  for  larger  separation  are 
also  consequences  of  the  charge-transfer  terms,  which  continue  to  provide  a  two-fold  redun¬ 
dant  explicitly  antisymmetric  representation  in  this  limit,  whereas  the  bands  of  eigenvalues 
approaching  unity  at  large  separation  arise  from  the  covalent  structures  in  the  representa¬ 
tion.  The  results  of  Figure  3  suggest  the  eigenvalues  of  the  metric  matrix  can  distinguish 
among  the  types  of  structures  present  in  the  representation  employed,  largely  confirm  the 
redundancy  of  covalent  and  ionic  terms  commonly  employed  in  molecular  calculations,  and 
emphasize  the  importance  of  removing  these  redundancies  in  large-basis-set  calculations. 

Illustrations  of  the  capabilities  of  the  Slater-based  spectral-product  codes  devised  for  energy 
calculations  are  given  in  Figures  4  and  5,  which  depict  low-lying  (n=l  and  2)  H2  singlet  and 
triplet  potential  energy  curves,  respectively,  obtained  from  full  configuration-interaction 
calculations  in  an  optimal  5s3p2dlf  valence  Slater  basis.27  These  results  are  in  general  but 
not  precise  accord  with  the  most  accurate  previous  calculations  available  for  the  indicated 
states.28  Of  particular  interest  are  the  shapes  of  the  B1!]^  and  E,F1E+  curves,  the  presence 
of  a  double  well  in  the  latter  curve  and  in  the  G,K1E+  curve,  and  the  large  barrier  in 
the  f  3nu  state.  The  long-range  nature  of  the  B 1 E+  and  E,F1E+  curves,  which  provide 

This  page  is  Distribution  A:  approved  for  public  release;  distribution  unlimited.  Page  14 


good  approximations  to  the  accurate  values,  is  commonly  but  subjectively  attributed  to 
contributions  from  charge-transfer  terms.  The  calculated  G.K 1  E+  and  f  1  IT v/  curves  provide 
reasonably  quantitative  descriptions  of  the  known  detailed  shapes  of  these  two  curves, 
but  are  not  in  precise  agreement  with  the  results  of  highly  accurate  calculations.28  Of 
course,  highly  accurate  calculation  of  the  excited-state  potential  energy  curves  of  H2  are 
generally  performed  employing  special  basis  sets  on  a  state-by-state  basis.  It  is  therefore 
satisfying  that  the  results  of  Figures  4  and  5  are  obtained  from  single  diagonalizations  of 
matrices  constructed  in  Slater-based  spectral-product  representations,  which  methodology 
is  applicable  to  other  molecules  more  generally. 

B.  Diatomic  CH  Molecule. 

The  foregoing  spectral-product  representations  of  molecules  containing  hydrogen  atoms 
(H2,  h3)  are  closely  related  to  corresponding  valence-bond  descriptions,  the  two  ap¬ 
proaches  differing  largely  in  the  normalization  conventions  employed  for  orbital-product 
states,  and  a  pre-diagonalization  of  the  atomic  hydrogen  Hamiltonian  in  the  spectral- 
product  development.  The  situation  is  quite  different  for  molecules  containing  many- 
electron  atoms.  The  differences  between  metric  or  Hamiltonian  matrices  constructed  in 
valence-bond  and  spectral-product  representation  can  be  illustrated  in  these  cases  employ¬ 
ing  the  well-studied  diatomic  CH  molecule  as  a  simple  example.  In  this  example,  both 
valence-bond  and  spectral-product  metric  and  Hamiltonian  matrices  are  constructed  em¬ 
ploying  identical  hydrogen  atom  Is  spin  orbitals  and  ls2(2s22p2  +  2s2p3  +  2 p4)  carbon 
atom  valence-shell  configurations,  the  latter  giving  rise  to  fifty-five  distinct  atomic  mul- 
tiplet  states.  The  valence-bond  states  made  from  these  orbital  configurations  are  in  the 
forms  of  standard  tableau  functions  which  are  employed  directly  in  constructing  the  re¬ 
quired  matrices,23  whereas  in  the  spectral-product  representation  the  hydrogen  atom  2Se 
and  carbon  atom  [(2)3Pe,  (2)1De,  (2)1Se,  5S°,  3D°,  3P°,  ^^D0,  3S°,  1P°]  multiplet  states  are 
constructed  and  the  spectral-product  matrices  of  Eqs.  (11)  to  (13)  are  evaluated  employ¬ 
ing  algorithms  and  codes  devised  explicitly  for  this  purpose.  Specifically,  the  expression 
:  R )  =  <&^bH\r  :  R )  •  V^^(il)  is  employed  to  affect  the  basis  transformation, 
where  ViY  ,r>(R)  is  a  non- unitary  matrix  constructed  from  an  ordered  product  of  corre- 
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sponding  atomic  H  and  C  transformation  matrices.  This  vb  — >  sp  transformation  approach 
is  employed  in  constructions  of  the  metric  and  Hamiltonian  matrices  in  the  spectral-product 
representation  from  their  forms  calculated  in  the  valence-bond  representation. 

In  Figure  6  are  shown  eigenvalues  of  the  doublet-  and  quartet-state  metric  matrix  for 
diatomic  CH  constructed  in  the  aforementioned  atomic  valence  configurations  using  an 
optimized  valence  orbital  basis.27  The  solid  lines  refer  to  eigenvalues  of  the  matrix  in  the 
orthonormal  spectral-product  representation  of  Eq.  (3),  whereas  the  dashed  lines  are  those 
obtained  in  the  same  configurational  basis  employing  the  standard  tableau  functions  of 
valence-bond  theory.23  The  two  sets  of  eigenvalues  of  Figure  6  are  evidently  quite  different, 
those  in  the  valence-bond  description  grouping  in  accordance  with  the  normalizations  of 
standard  tableau  functions  transforming  under  particular  irreducible  representations  of 
Snt  and  approaching  different  limits  at  large  atomic  separation.  By  contrast,  there  are 
only  three  distinct  highly  degenerate  eigenvalues  of  the  spectral-product  metric  matrix 
in  Figure  6,  which  are  seen  to  approach  unity  in  accordance  with  the  orthonormality 
of  the  spectral-product  representation.  The  small  number  of  distinct  eigenvalues  of  the 
latter  representation  relates  to  use  of  only  n  =  1  shell  hydrogen  and  n=2  shell  carbon 
multiplet  configurations  in  the  calculations  and  to  the  orthogonality  of  the  representation. 
Relatedly,  although  the  allowable  eigenvalue  interval  for  CH  is  formally  Q=7\ / (6!  1 !  1 ! ) — T, 
the  representation  employed  in  Figure  6  is  too  small  to  exhibit  the  full  redundancy  of  the 
basis  of  Eq.  (4)  in  this  case. 

Just  as  the  metric  matrices  obtained  in  the  valence-bond  and  spectral-product  representa¬ 
tions  generally  differ,  so  also  do  the  corresponding  Hamiltonian  matrices.  Nevertheless,  the 
energy  eigenvalues  obtained  from  the  two  representations  employing  identical  atomic  con¬ 
figurations  must  be  identical,  providing  a  useful  test  of  the  algorithms  and  codes  devised  to 
perform  the  present  calculations.  In  Figure  7  are  shown  low-lying  potential  energy  curves 
for  CH  obtained  in  (ls)2Se  atomic  hydrogen  and  (ls22s22p2)[3Pe,1De,1Se]  atomic  carbon 
representations.  The  two  sets  of  potential  energy  curves  obtained  from  the  valence-bond 
and  spectral-product  representations  are  found  to  be  identical,  but  to  provide  very  poor 
representations  of  the  known  accurate  values.13  Specifically,  although  the  ground  state  is 
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correctly  predicted  to  be  of  2 II  symmetry,  the  calculated  chemical  binding  energy  is  much 
too  small,  the  predicted  A2A  state  falls  above  the  B2E_  state,  and  no  other  chemically 
bound  states  are  predicted  by  the  representation. 

In  Figure  8  are  shown  the  results  of  calculations  similar  to  those  of  Figure  7,  now  including 
the  additional  atomic  carbon  n— 2  shell  multiplet  configurations  ls22s2p3  and  ls22p4 .  The 
former  configuration  gives  rise  to  the  important  5S°  and  3D°  atomic  carbon  states  which 
are  seen  to  have  significant  effect  on  the  calculated  potential  energy  curves.  In  particular, 
the  a4E_,  A2  A,  B2E_,  and  C2E+  state  curves  are  significantly  lower  and  in  the  accepted 
order,  a  bound  a4E_  curve  is  now  obtained,  and  the  ground  X2II  state  is  also  more  strongly 
bound.  The  dramatic  lowering  of  the  a  1 E “  state  curve  is  due  to  the  additional  4E_  con¬ 
figurational  state  function  arising  from  the  atomic  carbon  ls22s2p3  valence  configuration, 
which  is  spectrally  concentrated  in  and  dissociates  to  the  5S°  atomic  carbon  state,  whereas 
the  significant  lowering  of  the  A2A  and  B2E“  state  curves  correlating  with  the  excited 
xDe  atomic  carbon  state  is  due  to  contributions  in  the  spectral-product  representation 
from  the  high-lying  (ls22s2p3)3D°  atomic  carbon  state.  Although  population  analysis  can 
also  be  employed  to  identify  the  presence  of  important  configurations  in  calculated  wave 
functions,  the  atomic  spectral  concentration  of  particular  atomic  components  in  the  rep¬ 
resentation  is  particularly  meaningful  in  that  these  states  have  physical  significance,  and 
are  not  subjectively  based  on  particular  choices  of  basis  sets,  methods  of  calculation,  or 
analysis  techniques.  Moreover,  the  contributions  from  particular  atomic-state  products 
are  read  off  directly  from  the  eigenvectors  obtained  in  the  representational  basis  without 
further  manipulations. 

C.  Triatomic  CH2  Molecule. 

The  methylene  molecule  has  played  an  important  role  in  the  development  of  accurate 
quantum  chemical  methodologies,14  and  is  relevant  in  connection  with  the  metric  matrices 
and  atomic  spectral  compositions  of  molecular  wave  functions  of  interest  here.  Calculations 
similar  to  those  reported  above  for  the  CH  molecule  help  to  further  clarify  aspects  of  these 
issues,  employing  atomic  multiplet  configurational  models  for  illustrative  purposes. 
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In  Figure  9  are  shown  metric  matrix  eigenvalues  in  both  valence-bond  and  spectral-product 
representations  employing  the  hydrogen  atom  Is  and  carbon  atom  Is2  (2s22p2 +2s2p3 +2p4) 
configurations  described  in  the  preceding  Section.  In  this  case,  the  redundancy  factor  for 
the  corresponding  basis  of  Eq.  (4)  is  Q  —  8! / (6!  1!  1!)  =  56,  although  the  closed  Is2 
carbon  and  Is  hydrogen  shells,  and  the  absence  of  higher  orbital  excitations,  limits  the 
calculated  eigenvalues  to  a  much  smaller  spectral  interval.  As  in  Figure  6  for  the  CH 
molecule,  the  eigenvalues  of  Figure  9  are  very  different  in  the  two  representations,  with  the 
degenerate  spectral-product  values  uniformly  approaching  unity  at  large  atomic  separation 
and  those  in  the  valence-bond  representation  forming  distinct  groups  associated  with  the 
normalizations  of  standard  tableau  functions  transforming  under  the  different  irreducible 
representations  of  Sn± .  Of  course,  the  corresponding  Hamiltonian  matrices  in  the  two 
representations  also  differ,  although  the  energy  eigenvalues  must  be  identical  when  the 
same  orbital  configurations  are  employed  in  the  two  developments. 

In  Figures  10  and  11  are  shown  triplet-  and  singlet-state  potential  energy  curves  for 
methylene  in  symmetric  collinear  (H-C-H)  arrangements  obtained  in  the  multiplet  car¬ 
bon  (ls22s22p2)  and  hydrogen  (Is)  atomic  configurations  employed  in  the  foregoing.  The 
lowest-lying  potential  curves  in  this  minimal  valence  multiplet  representation  are  seen  to 
be  very  weakly  bound  3ng  and  1  Yl;j  states,  whereas  the  correct  lowest-energy  triplet  and 
singlet  states  in  symmetric  collinear  CH2  arrangement  are  known  to  be  the  X3£“  and  a1  Ag 
states.14  Additionally,  all  other  curves  in  Figures  10  and  11  are  seen  to  be  strongly  non¬ 
bonding.  Clearly,  the  minimal  multiplet  representation  fails  to  provide  even  qualitatively 
correct  descriptions  of  the  potential  energy  curves  in  CH2. 

In  Figures  12  and  13  are  shown  potential  energy  curves  for  the  states  of  Figures  10  and  11, 
now  constructed  employing  the  full  carbon  atom  ls2(2s22p2  +  2s2p3  +  2 p4)  valence-shell 
configurations  also  employed  in  the  preceding  CH  molecule  calculations.  As  in  the  case  of 
CH,  there  are  dramatic  changes  in  the  CH2  potential  curves  of  Figures  12  and  13  relative  to 
those  of  Figures  10  and  11.  Specifically,  the  X3£“  state  is  seen  to  become  the  lowest  energy 
state,  the  axAff  states  becomes  the  lowest-lying  single  state,  and  the  cdS4-  state  becomes 
the  first  excited  singlet,  in  accord  with  the  results  of  accurate  calculations. 14  The  lowering 
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of  the  X3X“  potential  curve  is  a  consequence  of  the  additional  3S~  configuration  arising 
from  the  presence  of  the  5S°  carbon  atom  state  in  the  larger  representation,  whereas  the 
lowering  of  the  a1Aff  and  cdS3"  states  is  due  to  molecular  states  associated  with  the  high- 
lying  3D°  carbon  atom  state  in  the  larger  atomic  representation  employed.  These  assertions 
are  based  on  examinations  of  the  eigenvectors  of  the  corresponding  states,  emphasizing  the 
convenience  of  the  atomic  spectral-product  representation  for  chemical  diagnostic  purposes. 

IV.  Discussion  and  Concluding  Remarks 

The  spectral-product  approach  to  ab  initio  molecular  electronic  structure  calculations 
would  seem  to  provide  a  fresh  perspective  of  sufficient  promise  to  warrant  further  effort  in 
its  computational  implementation.  The  theoretical  development  presented  in  Section  II  is 
largely  pedestrian  in  nature,  but  demonstrates  rigorously  the  formal  equivalence  between 
enforcement  of  electron  antisymmetry  prior  to  or  subsequent  to  Hamiltonian  matrix  eval¬ 
uation  in  the  atomic  spectral-product  representation.  In  the  case  of  prior  enforcement,  the 
totally  antisymmetric  but  redundant  representation  employed  gives  rise  to  Hamiltonian 
matrix  elements  each  of  which  formally  depends  on  the  positions  of  all  the  atoms  in  the 
aggregate,  requiring  repeated  evaluations  of  such  matrix  elements  in  construction  of  poten¬ 
tial  energy  surfaces.  Moreover,  the  linearly-independent  subspace  of  this  redundant  repre¬ 
sentation  must  be  isolated  to  avoid  encountering  computational  instabilities  in  obtaining 
molecular  energy  eigenfunctions  and  eigenvalues.  By  contrast,  the  post-ant isymmerization 
approach  employs  a  complete  but  not  over-complete  representation  which  gives  rise  to 
Hamiltonian  matrix  elements  which  can  potentially  be  evaluated  once  and  for  all  and 
retained  for  repeated  polyatomic  applications.  There  remains  the  task  in  this  approach 
of  isolating  the  totally  antisymmeric  subspace  of  the  spectral-product  representation  to 
obtain  the  physically  significant  molecular  eigenstates. 

The  matrix  representing  the  total  antisymmetrizer  in  the  spectral-product  basis,  which 
is  seen  from  the  development  of  Section  II  to  be  equivalent  to  the  metric  matrix  of  its 
antisymmetrized  form,  provides  a  method  for  isolating  the  totally  antisymmetric  subspace 
of  the  spectral-product  representation,  as  well  as  the  linearly-independent  subspace  of  its 
antisymmetric  form,  by  unitary  transformation  of  the  Hamiltonian  matrix.  This  unitary- 
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transformation  approach  to  antisymmetrization  provides  an  alternative  to  prior  global 
antisymmetrization  which  need  not  entail  direct  diagonalization  of  the  complete  aggregate 
metric  matrix.  Rather,  alternatives  can  be  devised  to  distinguish  among  atomic  interac¬ 
tions  for  which  antisymmetry  is  strictly  required  and  those  interactions  for  which  it  is  not 
required.  Since  only  overlap  integrals  are  needed  in  evaluations  of  the  metric  matrix,  with 
one-  and  two-electron  energy  integrals  not  needed  for  this  purpose,  there  is  considerable 
opportunity  to  develop  approximation  methods  which  can  take  advantage  of  the  atomic 
pairwise-interaction  nature  of  the  aggregate  Hamiltonian  matrix. 

The  computational  applications  to  simple  molecules  (H2,  H3,  CH,  CH2)  reported  in  Section 
III  illustrate  the  role  of  the  eigenvalue  spectrum  of  the  metric  matrix  in  simple  diatomic 
and  polyatomic  cases  in  assessing  closure  in  the  spectral-product  representation,  and  in  cor¬ 
respondingly  identifying  redundancy  in  its  explicitly  antisymmetrized  form.  The  allowable 
range  of  the  eigenvalue  spectrum  in  the  spectral-product  basis  [0,  Q  =  nt\/(ni\ri2\  ■  ■  -njv!)] 
is  seen  to  be  determined  by  the  number  of  electrons  in  each  atom  and  the  total  sum  of 
these  in  the  particular  normalization  of  the  antisymmetrizer  employed.  Convergence  of  the 
eigenvalues  to  the  extreme  points  of  the  allowable  range  (0  and  Q),  in  accordance  with  the 
compact-operator  nature  of  the  total  antisymmetrizer,  corresponds  to  the  separation  of 
totally  antisymmetric  and  non-totally  antisymmetric  states  in  the  spectral-product  basis, 
or,  equivalently,  to  the  separation  of  linearly-independent  and  linearly-dependent  states 
in  its  antisymmetrized  form.  The  introduction  of  additional  terms  in  the  representation, 
particularly  of  the  charge-transfer  terms  commonly  employed  in  molecular  calculations,  is 
seen  to  disrupt  this  simple  picture,  and  to  extend  the  allowable  range  of  the  eigenvalue 
spectrum  to  accommodate  the  additional  basis-set  redundancy  introduced.  It  seems  clear 
that  the  metric  matrix  provides  a  device  for  developing  suitable  atomic  spectral-product 
representations  entirely  in  the  absence  of  the  molecular  electronic  energy  calculations  com¬ 
monly  employed  in  devising  orbital  basis  sets  for  use  in  more  conventional  calculations. 

The  illustrative  molecular  energy  calculations  reported  in  Section  III  do  not  provide  defini¬ 
tive  quantitative  results  for  the  simple  molecules  studied  (H2,  CH,  CH2),  but  rather  demon¬ 
strate  to  some  degree  the  capabilities  of  the  valence-bond-  and  Slater-orbital-based  algo- 
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rithms  and  codes  devised  for  this  purpose,  as  well  as  the  nature  of  the  descriptor  of  molec¬ 
ular  electronic  structures  the  development  provides.  Calculations  of  the  twelve  excited- 
state  potential  energy  curves  in  H2  dissociating  to  the  n  —  2  atomic  limit,  obtained  from 
straightforward  diagonalizations  of  singlet  and  triplet  full  configuration-interaction  energy 
matrices  constructed  with  standard  Slater-orbital  basis  sets,  are  found  to  be  in  good  accord 
with  the  most  accurate  calculations  available,  as  obtained  largely  but  not  entirely  from  a 
series  of  individual  calculations  performed  with  special  basis  sets  not  applicable  to  other 
molecules  more  generally.  Additionally,  the  important  role  of  the  2s2p3  configuration  in 
atomic  carbon  is  understood  from  the  new  atomic  spectral-composition  perspective  pro¬ 
vided  by  the  spectral-product  development.  Specifcally,  the  contributions  to  molecular  CH 
and  CH2  eigenstates  from  5S°  and  3D°  atomic  carbon  states  which  arise  specifically  from 
the  2s2p3  configuration  are  seen  to  have  significant  qualitative  and  quantitative  effects  on 
the  positions  and  shapes  of  selected  potential  energy  curves.  Whereas  descriptions  of  such 
effects  in  terms  of  orbital  configurations  are  clearly  subjective,  the  perspective  provided  by 
the  spectral  compositions  of  molecular  electronic  wave  functions  in  the  spectral-product 
representation  is  not,  and  involves  the  physically  significant  many-electron  atomic  eigen¬ 
states  of  the  bonding  atoms.  It  would  seem  the  spectral-product  representation  provides 
both  a  conceptual  and  a  computational  basis  for  studies  of  molecular  electronic  eigenstates. 

Additional  calculations  not  reported  here  involving  H,  C,  N,  and  O  atom-containing  com¬ 
pounds  largely  support  and  amplify  the  foregoing  general  conclusions.  Studies  in  progress 
are  now  focused  on  development  of  signifcantly  more  efficient  and  robust  computational 
methods  for  performing  spectral-product  calculations,  use  of  larger  Slater-orbital  repre¬ 
sentations,  and  novel  means  of  isolating  the  required  invariant  totally  antisymmetric  sub¬ 
spaces.  Progress  in  these  areas  will  be  reported  elsewhere  in  due  course. 
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Figure  1.  Eigenvalues  of  the  doublet-state  metric  matrix  of  Eq.  (13)  for  symmetric 
(. Rab  —  Rbc )  collinear  H3  as  functions  of  atomic  separation  (Ha-H&-Hc):  Panel  (a)  2slp 
atomic  basis  (375  terms);  Panel  (b)  3s2p  atomic  basis  (2187  terms).  The  limiting  points 
of  the  spectra  are  0  and  Q  —  3!/ ( 1!  1!  1!)  =  6,  specifying  the  redundancy  of  the  basis  of  Eq. 
(4)  in  this  case. 
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Figure  2.  Const  ant- value  contours  for  the  two  doublet  eigenvalue  surfaces  of  the  metric 
matrix  of  Eq.  (13)  for  T-shaped  arrangements  in  H3  constructed  in  a  minimal  Is3  rep¬ 
resentation.  Panel  (a)  -  solid  lines  give  the  upper  surface;  Panel  (b)  -  solid  lines  give  the 
lower  surface,  employing  increments/decrements  of  0.1  between  adjacent  contours  in  each 
case.  The  cusps  in  the  contours  of  both  panels  indicate  the  presence  of  the  seam  of  crossing 
of  the  two  surfaces  along  the  w  —  1.155  h  line,  corresponding  to  D:>/,  symmetry.25 
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Figure  3.  Eigenvalues  of  the  singlet-  and  triplet-state  metric  matrix  of  Eq.  (13)  for 
diatomic  hydrogen  constructed  in  a  7s5p3d2flg  Sturmian  Slater  basis  set,2'  including 
both  covalent  (119)  and  ionic  or  charge-transfer  (154)  terms.  The  origins  of  the  structures 
evident  in  the  plot,  and  of  the  spectral  interval  (0,4)  in  this  case,  are  discussed  in  the  text. 
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Figure  4.  Singlet-state  potential  energy  curves  dissociating  to  n=l  and  2  limits  in  di¬ 
atomic  hydrogen,  obtained  from  full  (1953  terms)  configuration-interaction  calculations  in 
an  optimized  5s3p2cllf  valence  basis  of  Slater  orbitals.27 
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Figure  5.  Triplet-state  potential  energy  curves  dissociating  to  n=l  and  2  limits  in  di¬ 
atomic  hydrogen,  obtained  from  full  (1891  terms)  configuration-interaction  calculations  in 
an  optimized  5s3p2cllf  valence  basis  of  Slater  orbitals.27 
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Figure  6.  Eigenvalues  of  the  doublet-  and  quartet-state  metric  matrix  of  Eq.  (13) 
for  diatomic  CH,  constructed  employing  hydrogen  (Is)  and  carbon  ls2(2s22p2  +  2s2p3  + 
2 p4)  multiplet  atomic  configurations  in  an  optimized  valence  orbital  basis.27  The  solid 
lines  refer  to  eigenvalues  of  the  matrix  in  the  orthonormal  spectral-product  representation, 
whereas  the  dashed  lines  are  those  obtained  in  the  same  configurational  basis  employing 
the  normalization  conventions  of  the  standard  tableau  functions  of  valence-bond  theory.23 
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Figure  7.  Potential  energy  curves  for  diatomic  CH  constructed  employing  minimal 
valence-shell  multiplet  carbon  (ls22s22p2)  and  hydrogen  (Is)  atomic  configurations.  The 
spectroscopic  state  labels  employed  follow  conventions  from  experimental  and  previously 
reported  theoretical  studies.13 
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Figure  8.  Potential  energy  curves  for  diatomic  CH  obtained  in  valence-shell  multiplet 
carbon  ls2(2s22p2  +  2s2p 3  +  2 pA)  and  hydrogen  (Is)  atomic  configurations.  The  state 
designations  of  the  lower-lying  curves  are  as  in  Figure  7,  whereas  the  additional  curves 
dissociating  to  the  5S°  and  3D°  carbon-atom  limits  are  labeled  employing  conventional 
spectroscopic  notation. 13 
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Figure  9.  Eigenvalues  of  the  singlet-  and  triplet-state  metric  matrix  of  Eq.  (13)  for  sym¬ 
metric  collinear  (H-C-H)  methylene  arrangements,  constructed  employing  valence-shell 
multiplet  carbon  ls2(2s2p2  +  2s2p3  +  2 p4)  and  hydrogen  (Is)  atomic  configurations.  The 
solid  lines  refer  to  eigenvalues  of  the  matrix  constructed  in  the  orthonormal  spectral- 
product.  representation,  whereas  the  dashed  lines  are  those  obtained  in  the  same  configu¬ 
rational  basis  employing  the  normalization  conventions  of  the  standard  tableau  functions 
of  valence-bond  theory.23 
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Figure  10.  Potential  energy  curves  for  triplet  methylene  states  in  symmetric  collinear 
(H-C-H)  arrangements  obtained  in  minimal  valence-shell  multiplet  carbon  ( ls22s22p2)  and 
hydrogen  (Is)  atomic  configurations.  The  ground  X3S“  state  is  labeled  on  basis  of  exper¬ 
imental  and  previously  reported  theoretical  studies,14  whereas  the  other  states  are  given 
sequential  numerical  labels. 
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Figure  11.  Potential  energy  curves  for  singlet  methylene  states  in  symmetric  collinear 
(H-C-H)  arrangement  obtained  in  minimal  valence- shell  multiplet  carbon  ( ls22s22p2)  and 
hydrogen  (Is)  atomic  configurations.  The  two  states  given  spectroscopic  state  labels  follow 
conventions  from  experimental  and  previously  reported  theoretical  studies,14  whereas  the 
other  states  are  given  sequential  numerical  labels. 
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Figure  12.  Potential  energy  curves  for  triplet  methylene  states  in  symmetric  collinear 
(H-C-H)  arrangements  obtained  in  valence-shell  multiplet  carbon  Is2 (2s22p2  +  2s2p3  +  2p4) 
and  hydrogen  (Is)  atomic  configurations.  The  state  designations  of  the  curves  are  as  in 
Figure  10. 
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Figure  13.  Potential  energy  curves  for  singlet  methylene  states  in  symmetric  collinear 
(H-C-H)  arrangement  obtained  in  valence-shell  multiplet  carbon  ls2(2s22p2  +  2s2p3  +  2p4) 
and  hydrogen  (Is)  atomic  states.  The  state  designations  of  the  curves  are  as  in  Figure  11. 
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